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Precede de determination d'un bruit fluoroscopique 

L'invention a pour objet un precede de determination d'un bruit 
fluoroscopique. Le domaine de ('invention est celui de I'imagerie et plus 
precisement celui de la reduction du bruit fluoroscopique dans les images. Le 
domaine de l'invention est plus particulierement celui de la reduction du bruit 
dans les images acquises en sequences temporelies. 

Un but de l'invention est la reduction du bruit fluoroscopique dans des 

images. 

Un autre but de l'invention est la reduction du bmit fluoroscopique 
dans des images numeriques issues d'un precede d'acquisition 
fluoroscopique. 

Un autre but de {'invention est la determination d'une determination 
robuste d'un bruit fluoroscopique. 

Un autre but de l'invention est la determination d'une determination 
temps reelle d'un bruit fluoroscopique. 

Un autre but de l'invention est de reguler la quantite de rayonnement 
re9ue par un patient au cours d'un examen radiographique. 

Dans I'etat de la technique on sait qu'en fluoroscopie il existe un bruit 
fluoroscopique et que I'ecart type de ce bruit fluoroscopique est en racine 
carree du nombre de photons arrivant sur le detecteur. 

Dans I'etat de la technique, ia presence de ce bruit est connue mais 
n'est pas specifiquement traite autrement que par I'experience d'un 
observateur utiiisant les images pour en effectuer une interpretation. II se 
pose des lors un probleme d'interpretation, mais aussi un probleme 
d'exposition du patient. On peut en effet etre amene a augmenter la dose de 
rayonnement pour am6liorer le rapport signal a bruit de I'image. c'est a dire 
pour augmenter la proportion d'information pertinente dans I'image. Cette 
augmentation du rayonnement n'est pas sans consequence pour les cellules 
du patient qui risquent alors d'etre d^truites ou endommagees. 

Dans ['invention on resout ces problemes en filtrant i'image obtenue 
de maniere a ameliorer sa qualite. Ainsi, dans l'invention on elimine du bruit 
fluoroscopique apres I'avoir determine / modelise. Ce bruit est totalement 
determine, done modelise, par son ecart type. Get ecart type est lui-meme 
fonction de la racine carre du nombre de photon regu par le detecteur. Le 
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nombre de photon est lui lie au niveau de gris dans une image numerique. 
On utilise done des images numeriques pour realiser la modelisation. Cette 
modelisation est realisee en plusieurs etapes. Une premiere etapes consiste 
en I'acquisition de deux images numeriques d'une meme zone. Pour la suite 
5 une image est, sauf indication contraire, numerique. On note ici que 
I'enseignement de I'invention est valable quelque soit la nature, numerique 
ou analogique, du capteur permettant d'obtenir les images. Chaque pixel, ou 
point, d'une image est appaire avec un pixel de I'autre image via ses 
coordonnees dans Pimage. Chaque pixel a aussi une valeur de niveau de 

10 gris, ou niveau de gris. On regroupe les pixels par intervalle de niveaux de 
gris et on obtient ainsi des sous groupes de pixels. Pour chaque sous groupe 
de pixels on effectue une discrimination de la maniere suivante: on calcul la 
moyenne p et Tecart type a de Pi(x,y) - Pi-1(x,y), ou Pi(x, y) est le niveau de 
gris du pixel de Timage i et de coordonnees (x, y). Puis dans un sous groupe, 

15 on ne garde que les P(x,y) tel que Pi(x,y) - Pi-1(x,y) < p + k. o. On repete 
cette discrimination de maniere iterative sur le resultat de la discrimination 
precedente. Cette discrimination permet de reduire fortement le bruit de 
bouge. Pour obtenir des sous groupes encore plus fiable on elimine ceux qui, 
a la fin de la discrimination, ne sont pas centres, c'est a dire ceux dont la 

20 moyenne est superieure a 1.5 fois I'ecart type. On est alors en connaissance 
d'une collection de paire (v, a), ou v est un niveau de gris. Ces paires 
permettent d'effectuer une regression aboutissant a des parametres a, p et y 
tel que a(v) = a.Vv + p.v + y* ou o(v) est une modelisation du bruit 
fluoroscopique. Cette regression est rendue robuste en Titerant apres avoir 

25 pondere les a des couples (v, a) tirant la courbe vers le haut de maniere a 
obtenir une courbe passant sous la majorite des points (v, o). 

Une fois le bruit modelise, on peut le supprimer de Timage et ainsi 
ameliorer la lisibilite de I'image et eventuellement diminuer la dose de 
rayonnement requise pour Tobtention d'une image lisible. 

30 ^invention a done pour objet un precede de determination d'un bruit 

fluoroscopique comportant les etapes suivantes: 

- on acquiert au moins deux images i-1 , et i, d'une meme zone, 

- on transcode les images acquises en images numeriques identifiable 
a des matrices de dimensions horizontale x verticale egale a N x M, chaque 

35 image numerique etant alors composee de N.M points, chaque point d'une 
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image i etant identifiable par ses coordonnees 0 < x < N, et 0 < y < M, on 
parle alors d'un point Pi(x,y), a chaque Pi(x,y) correspond une valeur v 
resultat de Tacquisition de I'inaage, la valeur v ayant une dynamique comprise 
entre Vmin et Vmax, 
5 procede caracterise en ce qu'il comporte les etapes suivantes: 

- on divise la dynamique Vmax - Vmin en sous groupes definis par 
une borne inferieure Bi et une borne superieure Bs, les sous groupes etant 
d'intersection nulle» la reunion des sous groupes couvrant la dynamique 
Vmax - Vmin, un point d'une Image i appartient alors a un sous groupe 

10 donne lorsque Bi <= Pi(x.y) < Bs, 

- on calcule, pour au moins un sous groupe SG, I'ecart type a des 
Pi(x,y)-Pi-1(x,y), 

- on discrimine les Pi(x,y) de SG en ne conservant que ceux tels que 
le critere C : 

1 5 Pi(x,y) - Pi-1 (x,y) < p(Pi(x,y) - Pi-1 (x.y)) + k.a, 

soit satisfait, on obtient ainsi un sous groupe SG', 

- on applique les memes traitements au sous groupe SG' qu'au sous, 
groupe SG et ainsi de suite par iteration jusqu'a obtenir un sous groupe SG" 
correspondant a un critere de fin d'iteration, 

20 - on effectue les traitements iteratifs sur tous les sous groupes definis 

dans la dynamique Vmax - Vmin, pour chaque sous groupe on obtient ainsi 
un ecart type associe a une abscisse v = (Bi + Bs) / 2, 

- on effectue une regression sur les points obtenus a Tetape 
precedente pour determiner les coefficients a, p et y de la fonction de bruit: 

25 a(v) = a.Vv + (3.v + y 

definissant le bruit fluoroscopique pour un v donne. 

L'invention sera mieux comprise a la lecture de la description qui suit 
et a I'examen des figures qui I'accompagnent. Celles cl sont presentees a 
titre indicatif et nullement limitatif de Tinvention. Les figures montrent: 

30 - Figure 1: une illustration d'etapes du procede selon l'invention. 

- Figure 2: une illustration d'une sequence d'image. 

- Figure 3: une illustration d'une representation graphique des paires 

(V. a). 

La figure 1 montre une etape 101 preliminaire d'acquisition d'images. 
35 Cette acquisition est realisee par un appareil radiographique de type scanner 
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permettant d'acquerir des images d'organes internes d'un organisme vivant, 
et en particulier d'un corps humain. Les technologies de ces appareils sont 
basees sur remission d'un rayonnement qui est regu sur un capteur apres 
avoir traverse Torganisme. Le rayonnement regu par le capteur est done 
5 fonction de I'organisme traverse. Le capteur produit une image numerique 
susceptible d'etre visualisee sur un ecran, imprimee, et ou traitee par un 
ordinateur, ou tout appareil comportant des circuits de traitement. Le precede 
selon I'invention est enregistre dans une memoire sous forme de codes 
instructions et mis en oeuvre par un microprocesseur, ce en vue d'ameliorer 

10 rinterpretation d'une image et d'ou de diminuer la dose de rayonnement 
employe pour realiser un examen. Cette memoire et ce microprocesseur sont 
de preference compris dans Tappareil radiographique. Dans la pratique cette 
memoire et ce microprocesseur sont connectes a I'appareil radiographique, 
soit par un connecteur interne, soit par un connecteur exteme a Tappareil 

1 5 radiographique. 

Classiquement une image a une resolution de 1000 par 1000 et 
comporte done 1 millions de points ou pixels, chaque point de i'image etant 
repere par ses coordonnees dans I'image, chaque point de Timage etant 
associe a une valeur de niveau de gris, dintensite ou a une autre grandeur 

20 numerique. Dans notre exemple nous'considererons qu'il s'agit d'un niveau 
de gris obtenu apres exposition fluoroscopique, c'est a dire une exposition a 
un faible rayonnement. Une image est assimilable a un tableau de N lignes 
par M colonnes. N et M ne valent pas systematiquement 1000, leurs valeurs 
dependent de la resolution spatiale du capteur / detecteur utilise. Leurs 

25 valeurs n'infiuent pas sur Tefficacite du precede. Le mode operatoire de 
i'invention reste valable pour une exposition non fluoroscopique. 

Comme deja dit, en fluoroscopie I'ecart type du bruit fluoroscopique 
est en racine carree du nombre de photons arrivant sur le detecteur, hors il 
se trouve que le niveau de gris est proportionnel a la quantite de photon 

30 regu, ce qui permet de travailler sur les niveaux de gris. 

Le bruit fluoroscopique est la resuitante du bruit quantique et de la 
reponse du systeme d'acquisition d'image c'est a dire du detecteur, Le 
precede decrit par la suite traite I'ensemble bruit quantique et appareil. 
Techniquement parlant, le precede n'est pas specifique au bruit quantique ou 

35 a un appareil, mais il peut etre appliquer appliquer a tout autre bruit non 
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spatialement correle. 

Pour la plupart des examens les images sont acquises en sequences 
temporelles. c*est a dire les unes a la suite des autres. Cela permet 
d'observer un organe au cours d'un cycle comme ie cycle cardiaque par 
5 exemple. Dans de tel examen on dispose done de plusieurs images, ou 
cliches, successives d'une meme zone d'un organisme. Ces images 
successives ont la meme resolution et deux pixels de meme coordonnee 
correspondent a une meme zone de I'organisme, sauf si ('organisme s'esl 
deplace entre les deux cliches. De tels deplacements peuvent survenir 
10 quand on essaie d'obtenir une image du coeur par exemple. 

Considerons que Tappareil radiographique a fourni deux images i-1 et 
i d'une meme zone d*un organisme. Ces images ont une resolution horizontal 
de N pixels, et une resolution verticale de M pixels. Comme illustre par la 
figure 2. Un point de coordonnees (x, y) d'une image i a un niveau de gris 
15 egal a Pi(x, y). 

De Tetape 101 on passe a une etape 102 de sous-echantillonage des 
images issues de Tappareil de radiographie. A chaque pixel d'une image est 
associe une valeur correspondant a un niveau de gris. Cette vaieur a une 
certaine dynamique. En d'autre termes a chaque pixel est associee une 
20 valeur comprise entre Vmin et Vmax. Classiquement Vmin vaut 0 et Vmax 
2'^14. Dans la pratique Vmin et Vmax dependent du capteur utilise et de son 
calibrage. Lorsque Ton sous echantillonne une image, on divise Tintervalle 
[Vmin, Vmax] en sous intervalles disjoints de telle maniere que la reunion 
des sous intervalles couvre [Vmin, Vmax]. A chaque sous intervalle [Bi, Bs] 
25 inclut dans [Vmin, Vmax] correspond un sous groupe de pixels tels que Bi <= 
Pi(x,y) < Bs. 

Dans la pratique les sous inten/alles ont une dynamique de 20, c'est a 
dire que Tintervalle [Vmin, Vmax] est decoupe en sous intervalle de largeur 
20. Le sous echantillonnage consiste alors a repartir les points d'une image 
30 dans les sous groupes correspondant aux sous intervalles. II est possible 
qu'il y ait des sous groupes vides. Un sous groupe est aussi appele un 
echantillon. 

Dans la pratique le sous echantillonnage est realise des que Ton a 
decide du decoupage en sous intervalles et il n'y a pas de repartition 
35 physique. La repartition se fait par simple lecture du niveau de gris d'un pixel. 
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La repartition en sous groupe se fait par lecture de llmage i. 

De I'etape 102 on passe a une etape 103 de calcul d'un ecart type 
pour un echantillon. On considere les points d'un sous groupe SG de Timage 
i en correspondance avec des points de I'image i-1. Des points sont en 
correspondance des lors qu'ils ont les menne coordonnees. Ainsi le fait 
d'avoir determine un sous groupe dans Timage i determine un sous groupe 
comportant les meme points dans Timage i-1. Dans I'etape 103 on calcule 
done recart type a des (Pi(x, y) - Pi-1(x, y)) avec (x, y) appartenant a SG. 

De I'etape 103 on passe a une etape 104 de discrimination du sous 
groupe SG. Cette discrimination consiste a eliminer les points aberrants 
resultant par exemple d'un bouge de Torgane a imager. La discrimination se 
fait selon le critere suivant: (x, y) de SG appartient a SG' si et seulement si 

Pi(x,y) - Pi-1(x,y) < |j(Pi(x,y) - Pi-1(x.y)) + k,o, 
ou M(Pi(x,y) - Pi-1(x,y)) est la moyenne des (Pi(x,y) - Pi-1(x,y)) pour (x, y) 
dans SG. SG' est done inclus dans SG mais peut ne pas comporter tous les 
points de SG. 

Le nombre k est un parametre de I'algorithme, il vaut de preference 2. 
Mais n'importe quel nombre non nul eonvient. Plus le nombre k choisi est 
eleve moins I'algorithme est discriminant. Un nombre k superieur a 10 
eonvient done aussi mais la discrimination est alors tres faible. 

De I'etape 104 on passe a I'etape 105 de determination de la fin de 
I'iteration, Un premier critere possible pour la fin de I'iteration est que SG' soit 
identique a SG. Dans ce cas on passe a une etape 106 de centrage. Sinon 
on remplace SG par SG' et on reprend les etapes 103 et 104. 

Un deuxieme critere possible de fin d'iteration est un nombre de 
boucle. Dans ce cas, dans I'etape 105 on remplace SG par SG' est on 
reprend a I'etape 103. Ce remplacement est realise un nombre maximal 
predefini de fois par exemple 5 fois. Mais ce nombre est indicatif et 
parametrable en fonction de la qualite de la modelisation que Ton souhaite 
obtenir. Plus ce nombre est eleve, meilleur est la qualite, mais plus longue 
est la modelisation. 

Dans la pratique il est possible de combiner les premier et deuxieme 
criteres via une combinaison logique. Le critere est alors que Ton itere tant 
que le nombre d'iteration est inferieur au nombre maximal ET tant que SG 
est different de SG'. 
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Dans Tetape 105, si le critere de fin d'iteration n'est pas valide, on 
passe a I'etape 103 apres avoir effectuer I'operation SG = SC. Sinon, on 
passe a I'etape 106. 

L'etape 106 a pour but de rendre le procede encore plus robuste. 
5 Dans I'etape 106 on elimine du traitement les echantillons, ou sous groupe 
SG', resultant de I'etape 105 tel que: 

M(Pi(x,y) - PHl(x,y)) > l.a(Pi(x.y) - Pi-1(x,y)), 
avec (x, y) dans SG' et 1 valant classiquement 1.5. Le fait d'eliminer un sous 
groupe signlfie que ce sous groupe ne sera plus pris en compte a partir de 
10 I'etape 108 incluse. 

De I'etape 106 on passe a une etape 107 dans laquelle on verifie que 
tous les sous groupes determines lors de I'etape 102 ont ete traites. Si ce 
n'est pas le cas on repete les etapes 103 a 106 pour tous les sous groupes 
non encore traites. 

15 A la fin de I'etape 107, on est capable d'associer a chaque sous 

groupe non vide et non ecarte a I'etape 106 une valeur v et un ecart type o. 
La valeur v est la valeur moyenne du sous intervalle ayant servi a determiner 
initialement le sous groupe. Un sous intervalle etant determine par une borne 
inferieur Bi et une borne superieur Bs, v vaut (Bi + Bs) / 2. L'ecart type o est 

20 le dernier ecart type calcule pour le sous groupe. On a done une collection 
de couple (v, a). 

De I'etape 107, lorsqu'il n'y a plus de sous groupe a traiter, on passe a 
une etape 108 de regression. La figure 3 illustre le fait qu'il est possible de 
representer graphiquement les couples (v, a), en plagant v en abscisse, v est 
25 compris entre Vmin et Vmax. 

A partir des couples (v, o) on effectue une regression de maniere a 
determiner des premiers coefficients a, p et y modelisant une courbe o(v) = 
a.Vv + p.v + Y passant au plus pres de Tensemble des couples (v, a). 

De I'etape 108 on passe a une etape 109 dans laquelle on evalue un 
30 critere de fin d'iteration. Un premier critere de fin d'iteration est, par exemple, 
la presence de points sous la courbe definie a I'etape 108. Si de tels points 
existent on passe a une etape 110 dans laquelle on modifie la ponderation 
de ces points. De I'etape 110 on passe alors a I'etape 108. Si de tels points 
n'existent pas. on passe de I'etape 109 a une etape 1 1 1 de fin. L'etape 108 
35 est done effectuee plusieurs fois sur la base d'une collection de couple dont 
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la ponderation evolue en fonction de la position du couple par rapport a la 
derniere courbe determinee, c'est a dire par rapport aux derniers coefficients 
a, (3 et Y determines a I'etape 108. 

Un deuxieme critere de fin d'iteration est Tutilisation d'un compteur 
5 d'iteration. Avec ce deuxieme critere on effectue les etapes 108 a 110 un 
nombre R predetermine de fois. Dans la pratique R vaut 5. II est toutefois 
possible d'utiliser n'importe quelle valeur pour R. 

Dans une variante on utilise les premiers et deuxieme criteres 
combines par un operateur logique ce qui permet de prendre en compte la 
10 plus favorable des conditions d'arret de iteration en terme de temps de 
traitement. Ainsi on arrete I'iteration des que le nombre maximal d'iteration 
est atteint ou des qu'il n'y a plus de couple au dessous de la derniere courbe 
determinee. 

Dans I'etape 110 on prepare un afflnage de la regression de I'etape 

15 108. Cette affinage se fait via une ponderation des couples (v, a) situes au 
dessus de la courbe definie a I'etape 108. Ces points sont detectes 
facilement, il s'agit des points pour lesquels a(v) < a. Pour ces couples on 
applique une ponderation P strictement inferieure a 1 sur la valeur a du 
couple (v, a) qui devient alors (v, P.a) . De I'etape 1 10 on passe a I'etape 108 

20 de regression avec I'ensembles des couples (v, a), qu'ils soient ponderes ou 
non, issus de I'etape 110. Une valeur typique de P est 0,95, mais nimporte 
quelle valeur comprise dans Tintervalle [0 ... 0,99] convient. Des resultats 
satisfaisant sont obtenus avec une valeur P dans Tintervalle [0,75 ... 1]. 

Lorsque que la condition de fin d'iteration est remplie, on passe de 

25 I'etape 109 a I'etape 111 de fin. Dans I'etape 111 I'appareil ayant effectue le 
traitement est en possession de coefficients a, (3 et y de la fonction: 

a(v) = a.Vv + p.v + Y 
modelisant le bruit fluoroscopique lors d'un examen mettant en oeuvre un 
appareil fluoroscopique. L'appareil ayant effectue le traitement utilise ces 

30 coefficients soit pour determiner un filtre permettant de produire une image 
dans laquelle bruit fluoroscopique a ete sensiblement reduit, soit il utilise ces 
coefficients pour parametrer I'appareil fluoroscopique. 

Dans le cas d'un filtrage, I'image apres filtrage est done exempte de 
bruit fluoroscopique ce qui favorise sa lecture et son interpretation par un 

35 praticien. II s'agit done d'une avance considerable pour le praticien qui peut 
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alors se fier d'avantage aux images fournies par les appareils de 
radiographie. 

Dans le cas d*un parametrage de Tappareil fluoroscopique, cela est 
utile car une fois les coefficients determine, on est en connaissance de la 
5 quantite de bruit presente dans une image, et done du rapport signal sur 
bruit. Si ce rapport signal sur bruit est satisfaisant, c'est a dire beaucoup de 
signal et peut de bruit, on diminue la dose de rayonnement emise par 
I'appareil fluoroscopique. Uexamen devient alors moins traumatisant pour la 
patient. D'un autre cote, si le rapport signal sur bruit est mauvais on 
10 augmente la dose de rayonnement afin que Texamen soit pertinent. Dans les 
deux cas on obtient un gain car on a reussi a determiner la dose optimale de 
rayonnement pour obtenir le resultat le plus pertinent possible, 

Dans la pratique les deux utilisations des coefficients a, (3 et y sont 
mises en oeuvre. 

15 Une sequence d'image comporte un nombre variable d'images. On 

utilise les deux premieres images de la sequence pour determiner les 
coefficients a, (3 et y qui seront appliques par la suite pour le traitement de 
toute les images de la sequence. 

Dans une variante, on calcul les coefficients a, |3 et y pour chaque 

20 image de la sequence, chaque image est alors associee a son propre jeux 
de coefficients. 

Dans la description on utilise une image i pour definir les sous 
groupes, et on appaire Timage i avec Timage precedente i-1. Dans une 
variante on peut aussi appairer I'image i avec une image suivante 1 + 1. 

25 Le procede selon Tinvention est mis en oeuvre soit sur une station de 

traitement d'images numeriques, soit dans un dispositif de commande d'un 
appareil de radiographie. Avec les machines usuellement utilisees dans le 
domaine de la radiographie le procede selon invention fournit des resultats 
en une trentaine de millisecondes pour le traitement d'une image ayant une 

30 definition d'un million de points. Ces performances sont tres satisfaisantes 
dans un contexte temps reel. Un operateur n'a pas Timpression d'attendre 
I'image. 
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alors se fier d'avantage aux images fournies par les appareils de 
radiographie. 

Dans le cas d'un parametrage de I'appareil fluoroscopique, cela est 
utile car une fois les coefficients determine, on est en connaissance de la 
quantite de bruit presente dans une image, et done du rapport signal sur 
bruit. Si ce rapport signal sur bruit est satisfaisant. c'est a dire beaucoup de 
signal et peut de bruit, on diminue la dose de rayonnement 6mise par 
I'appareil fluoroscopique. L'examen devient alors moins traumatisant pour la 
patient. D'un autre cote, si le rapport signal sur bruit est mauvais on 
augmente la dose de rayonnement afin que l'examen soit pertinent. Dans les 
deux cas on obtlent un gain car on a reussi a determiner la dose optimale de 
rayonnement pour obtenir le resultat le plus pertinent possible. 

Dans la pratique les deux utilisations des coefficients a, (3 et y sont 
mises en oeuvre. 

Une sequence d'image comporte un nombre variable d'images. On 
utilise les deux premieres images de la sequence pour determiner les 
coefficients a, p et y qui seront appliques par la suite pour le traitement de 
toute les images de la sequence. 

Dans une variante, on calcul les coefficients a, (3 et y pour chaque 
image de la sequence, chaque Image est alors associee a son propre Jeux 
de coefficients. 

Dans ia description on utilise une image i pour definir les sous 
groupes, et on appaire I'image i avec I'image precedente i-1. Dans une 
variante on peut aussi appairer I'image i avec une image suivante 1+1. 

Le precede selon I'invention est mis en oeuvre soit sur une station de 
traitement d'images numeriques, soit dans un dispositif de commande d'un 
appareil de radiographie. Avec les machines usuellement utilisees dans le 
domaine de la radiographie le precede selon I'invention fournit des resuitats 
en une trentaine de millisecondes pour le traitement d'une image ayant une 
definition d'un million de points. Ces performances sont tres satisfaisantes 
dans un contexte temps reel. Un operateur n'a pas I'impression d'attendre 
I'image. 

Dans une variante le critere de fin d'iteration est un nombre 
d'iterations superieur a 5. 

Dans une variante R est dans I'intervalle [3 ... 10]. 
Dans une variante R est superieur ^10. 
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REVENDICATIONS 

1 - Procede de determination d'un bruit fluoroscopique comportant les 
etapes suivantes: 

5 - on acquiert (101) au moins deux images i-1, et i, d'une meme zone, 

- on transcode les images acquises en images numeriques identifiable 
a des matrices de dimensions horizontale x verticale egale a N x M, chaque 
image numerique etant alors composee de N.iVI points, chaque point d'une 
image i etant identifiable par ses coordonnees 0 < x < N, et 0 < y < M, on 

10 parle alors d'un point Pi(x,y), a chaque Pi(x,y) correspond une valeur v 
resultat de Tacquisition de Timage, la valeur v ayant une dynamique comprise 
entre Vmin et Vmax, 

procede caracterise en ce qu'il comporte les etapes suivantes: 

- on divise (102) la dynamique Vmax - Vmin en sous groupes definis 
15 par une borne inferieure Bi et une borne superieure Bs, les sous groupes 

etant d'intersection nulle, la reunion des sous groupes couvrant la dynamique 
Vmax - Vmin, un point d'une image i appartient alors a un sous groupe 
donne lorsque Bi <= Pi(x,y) < Bs, 

- on calcule (103) , pour au moins un sous groupe SG, I'ecart type a 
20 des Pi(x,y) - PI-1 (x.y), 

- on discrimine (104) les Pi(x,y) de SG en ne conservant que ceux tels 
que le critere C : 

Pi(x,y) - Pi-1(x,y) < |j(Pi(x,y) - Pi-1(x,y)) + k.a. 
soit satisfait, on obtient ainsi un sous groupe SG', 
25 - on applique (105) les memes traitements au sous groupe SG' qu'au 

sous groupe SG et ainsi de suite par iteration jusqu'a obtenir un sous groupe 
SG" correspondant a un critere de fin d'iteration, 

- on effectue (107) les traitements iteratifs sur tous les sous groupes 
definis dans la dynamique Vmax - Vmin, pour chaque sous groupe on obtient 

30 ainsi un ecart type associe a une abscisse v = (Bi + Bs) / 2, 

- on effectue (108) une regression sur les points obtenus a Tetape 
precedente pour determiner les coefficients a, p et y de la fonction de bruit: 

a(v) = a.Vv + p.v + Y 
definissant le bruit fluoroscopique pour un v donne. 
35 2 - Procede selon la revendication 1, caracterise en ce que avant la 
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regression, et apres la discrimination, on elimine (106) pour la suite du 
traitement les sous groupes non centres, c'est a dire les sous groupes tel 
que la moyenne du sous groupe est superieure a I fois I'ecart type, I valant 
de preference 1.5. 

5 3 - Procede selon Tune des revendications 1 ou 2 caracterise en ce 

qu*on applique la fonction de bruit aux images i pour reduire !e bruit 
fluoroscopique dans cette image. 

4 - Procede selon Tune des revendications 1 a 3, caracterise en ce 
que k est un nombre non nul. 
10 5 - Procede selon I'une des revendications 1 a 4, caracterise en ce 

que le critere de fin d'iteration est un nombre d'iterations superieur a 5, 

6 - Procede selon Tune des revendications 1 a 5, caracterise en ce 
que le critere de fin d'iteration est le fait que tous les points de SG' satisfont 
au critere C. 

15 7 - Procede selon I'une des revendications 1 a 6, caracterisee en ce 

que: 

- lors d'une premiere regression on obtient des premiers coefficients a, 
P et Y determinant une courbe separant les sous groupes en deux, ceux dont 
I'ecart type est au dessus de la premiere courbe, et ceux dont I'ecart type est 

20 en dessous de la premiere courbe, 

- on applique (109 - 110) une ponderation P inferieure a 1 a Tecart 
type des sous groupes dont I'ecart type est situe au dessus de la premiere 
courbe, 

- on effectue une deuxieme regression a partir des sous groupes 
25 ponderes pour obtenir des deuxiemes coefficients a', p' et y' determinant une 

nouvelle courbe de bruit, 

- on effectue a partir de la nouvelle courbe les meme traitement que 
ceux effectues a partir de la premiere courbe et ainsi de suite un nombre R 
de fois. 

30 8 - Procede selon la revendication 7, caracterise en ce que P est dans 

rintervalle [0.75 ... 0.99]. 

9 - Procede selon la revendication 7, caracterise en ce que P est dans 
rintervalle [0 ... 0.75]. 

10 - Procede selon I'une des revendications 7 a 9, caracterise en ce 
35 que R est dans rintervalle [3 ... 10]. 
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1 1 - Procede selon Tune des revendications 7 a 9. caracterise en ce 
que R est superieur a 10, 
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